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ABSTRACT 



Context. The progenitor mass of type IIP supernova can be determined from either hydrodynamic modeling of the event or pre-explosion 
observations. 

Aims. To compare these approaches, we determine parameters of the sub-luminous supernova 2005cs and estimate its progenitor mass. 
Methods. We compute the hydrodynamic models of the supernova to describe its light curves and expansion velocity data. 
Results. We estimate a presupernova mass of 17.3 ± 1 Mq, an explosion energy of (4. 1 ± 0.3) x 10^" erg, a presupernova radius of 600 ±140 Rq, 
and a radioactive '*Ni mass of 0.0082 ± 0.0016 Mq. The derived progenitor mass of SN 2005cs is 18.2 ± 1 Mq, which is in-between those of 
low-luminosity and normal type IIP supernovae. 

Conclusions. The obtained progenitor mass of SN 2005cs is higher than derived from pre-explosion images. The masses of four type IIP 
supernovae estimated by means of hydrodynamic modeling are systematically higher than the average progenitor mass for the 9-25 Mq mass 
range. This result, if confirmed for a larger sample, would imply that a serious revision of the present-day view on the progenitors of type IIP 
supernovae is required. 

Key words, stars: supernovae: individual: SN 2005cs - stars: supernovae: general 



1. Introduction 

Type IIP supernovae (SNe IIP) originate presumably from 
main-sequence stars of the mass range of 9 - 25 Mq (Heger 
et al. I2003D . If this is the case, the predicted rate of SNe IIP 
should follow the Salpeter initial mass function with half of 
events occurring for stars of mass below 13 Mq. This conjec- 
ture requires confirmation by means of the determination of 
progenitor mass for an extended sample of SNe IIP. 

There are two ways to recover the progenitor mass of 
SN IIP on the main sequence. The first method is detection 
of the presupernova (pre-SN) in archival images of the host 
galaxy. The estimated flux and color index of the detected pre- 
SN is then converted into a stellar mass using the flux and color 
index predicted by stellar evolution models. Data for the avail- 
able directly identified progenitors in the compilation of Li et 
al. (2007 ) indicate that for eight SNe IIP the progenitor masses 
have been estimated in this way, and for six SNe IIP, upper 
limits to the progenitor mass have been found. 

An alternative approach to the mass determination involves 
hydrodynamic modeling of light curves and expansion veloci- 
ties for the well-observed SNe IIP. Combining the ejecta mass 
derived from the hydrodynamic modeling with the mass of 
the neutron star and the mass lost by the stellar wind pro- 
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vides us with the mass estimate of the progenitor. Henceforth, 
the progenitor mass determined by this method is referred to 
as the "hydrodynamic mass". At present, the hydrodynamic 
mass is measured only for three SNe IIP: the peculiar type 
IIP SN 1987A (Woosley 1988; Blinnikov et al.l IOOOl Utrobin 
|2005] |. the normal type IIP SN 1999em (Baklanov et al. l2005l 
Utrobin |2007] |. and the low-luminosity type IIP SN 2003Z 
(Utrobin et al. l2007D . The small amount of SN IIP events with 
the measured hydrodynamic mass is related to the fact that the 
hydrodynamic modeling requires a complete multi-band pho- 
tometry at both the plateau and the radioactive tail, and spectra 
of sufficient quality. 

There are only a few other SNe IIP that also meet these 
requirements. Among these is the sub-luminous type IIP 
SN 2005cs (Pastorello et al. I2006I I. On the basis of its lumi- 
nosity and expansion velocities, this SN is intermediate be- 
tween the low-luminosity and normal SNe IIP. Parameters of 
SN 2005cs are of considerable interest for two major reasons: 
(1) this SN is expected to have intermediate parameters, which 
would be interesting to check; (2) the pre-SN was detected in 
the pre-explosion images and its progenitor mass was estimated 
by several groups (Maund et al. 120051 Li et al. 120061 Eldridge 
et al. 2007 ) providing an opportunity to compare different mass 
estimates. 
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In the paper, we perform hydrodynamic modeling of SN 
2005cs to recover the parameters: ejecta mass, explosion en- 
ergy, pre-SN radius, and radioactive ^^Ni mass. We start with 
the description of the model and observational data used and 
then present the results of the hydrodynamic modeling for 
SN 2005cs (Sect. |2]i. In Sect. [3] we present additional argu- 
ments in favor of our choice of pre-SN models. Possible un- 
certainties in the hydrodynamic mass of SN 2005cs are ana- 
lyzed (Sect. |4li, and finally the implications of hydrodynami- 
cally studied objects for the origin of SNe IIP are discussed 
(Sect.©. 

We adopt an explosion date on June 27.5 UT (JD 2453549) 
and a distance of 8.4 Mpc following Pastorello et al. (2006) . 
and a reddening E(B - V) = 0.12 taken from Li et al. ([2006] l. 

2. Optimal hydrodynamic model 

The hydrodynamic model applied to SN 2005cs is essentially 
the same as used before for SN 1999em (Utrobin 2007). The 
pre-SN structure is set to be a non-evolutionary model of a red 
supergiant (RSG) star (but see Sect. [3]). The chemical compo- 
sition of the hydrogen envelope is solar. Although this might 
be a simplification, there is no observational evidence that the 
hydrogen abundance in the atmosphere of the RSG, e.g. a Ori, 
differs notably from solar (Harper et al. l2001l l. The effect of the 
variation in the surface abundances on the SN IIP light curve 
was studied before, and it was found that even a significant 
change in the hydrogen abundance of outer layers only weakly 
affects the light curve (Utrobin 2007). The light curve of the 
SN 2005cs model of hydrogen abundance X=0.65 and helium 
abundance Y=0.33, predicted for an 18 Mq star with initial so- 
lar composition (Heger [T998l l. is almost indistinguishable from 
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Fig. 1. Density distribution as a function of interior mass a) and 
radius b) for the optimal pre-SN model. The central core of 1 .4 
Mq is omitted. 
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Fig. 2. The mass fraction of hydrogen (solid line), helium (long 
dashed line), heavy elements (short dashed line), and radioac- 
tive ^^Ni (dotted line) in the ejecta of the optimal model. 



the case of solar composition. The inner layers of the ejecta are 
mixed with the helium core in the same way as in the model for 
SN 1999em. 

We refer to the hydrodynamic model as being the optimal 
one, in terms of an "eye-fit" to the observational light curve and 
the evolution of photospheric velocity. In general, a numeri- 
cal optimization procedure could be developed to complete the 
search for the best-fit model. However, at present it would re- 
quire an enormous amount of computational time that would 
be unjustified because the error of eye-fit is far less than the er- 
ror introduced by uncertainties in the distance and interstellar 
extinction. The search for the optimal model of SN 2005cs uses 
the hydrodynamic model behavior in parameter space studied 
earlier in detail (Utrobin |2007l l. The ^^Ni mass is determined 
empirically from the comparison of the /?-band luminosity of 
SN 2005cs at the radioactive tail with that of SN 1987 A. For 
the adopted distance and reddening, the iJ-values at the age of 
250-300 days (Tsvetkov et al. |2006] | correspond to the ^''Ni 
mass of 0.0082 M0. 

The observed bolometric light curve of SN 2005cs is recov- 
ered from U BVRI photometry (Pastorello et al. 120061 Tsvetkov 
et al. .2006) using a black-body approximation for the SN radi- 
ation, while expansion velocities at the photosphere were taken 
from Pastorello et al. (120061 1. 

The dependence of the light curve at the initial adiabatic 
cooling stage and the plateau phase on SN parameters provides 
us with a toolkit to search for the optimal model in parameter 
space. The obtained model of SN 2005cs is characterized by 
the ejecta mass Mg„v, = 15.9 Mq, the explosion energy E - 
4.1 X 10^° erg, and the pre-SN radius - 600 Rq with the 
^^Ni mass M^i = 0.0082 Mq. The ejecta mass combined with 
the neutron star mass results in the pre-SN mass of 17.3 Mq. 

The optimal pre-SN density structure is shown in Fig. [T] 
and the chemical composition in Fig. |2l The helium core is 
mixed with the hydrogen envelope so that the hydrogen abun- 
dance increases linearly with mass in the inner 9 Mq. We adopt 
the helium-core mass of 5.4 Mq, which corresponds to the 
~ 18 Mq progenitor (Hirschi et al. 120041) . We note that the 
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basic SN parameters are insensitive to the helium-core mass 
(Utrobin et al. |2007l l. In the freely expanding envelope, the hy- 
drogen is mixed downward to 300 km s \ while the radioactive 
^^Ni is mixed outward to 610 kms"' (Fig.|3]l. 

The observed bolometric light curve is reproduced by our 
optimal model (Fig. |4|i. A small disparity between model and 
observations at the radioactive tail is probably caused by the 
black-body approximation applied to the reconstruction of the 
observed bolometric light curve. This approximation is cer- 
tainly rough at the nebular epoch. We note that the ^^Ni mass 
was derived from the comparison of the /?-band luminosities of 
SN 2005cs and SN 1987A at similar nebular epochs. 

Although the one-group approximation for radiation trans- 
fer prevents us from achieving a detailed description of the 
spectral energy distribution, our model does reproduce the gen- 
eral magnitude of the broad-band photometry, although the 
match is not precise. Both the observed and calculated B light 
curves (Fig. [5^) show an initial peak related to the cooling of 
hot outer layers of shocked ejecta. The amplitude and width of 
this peak is sensitive to the structure of the outermost rarefied 
layers of the pre-SN envelope, and the density structure shown 
in Fig.[T]is optimal in this sense. The later (f > 50 days) behav- 
ior of the flux in B band is poorly reproduced because the black- 
body approximation of the model spectrum in the A < 4500 A 
range is too crude at the plateau stage. Fortunately, this discrep- 
ancy does not affect the model fit to the bolometric light curve 
because of a small contribution of the B band to the bolomet- 
ric luminosity at this stage. The calculated V andR light curves 
for this model describe satisfactorily observations (Figs.lSji and 
c). The computed evolution of the photospheric velocity is also 
consistent with observations (Fig.|5}i), supporting the detailed 
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Fig. 3. The density and the ^^Ni mass fraction as a function of 
the velocity for the optimal model at f = 50 days. Dash-dotted 
line is the density distribution fit p oc v"^-^. 
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Fig. 4. Comparison of the calculated bolometric light curve 
of the optimal model {solid line) with the bolometric data of 
SN 2005cs evaluated from the photometric observations of 
Pastorello et al. ( |2006l l {open triangles) and Tsvetkov et al. 
(120061) {open circles). 

Table 1. Hydrodynamic models for evolutionary presuper- 
novae. 



Model 
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M pre-SN 


X Y 




(Rq) (Mo) (10^" erg) (IQ-^Mq) (Mq) 


(Mq) 




EMI 


1200 15.9 


4.1 


0.82 


5.4 


17.3 


0.65 0.33 


EM2 


600 15.9 


4.1 


0.82 


5.4 


17.3 


0.65 0.33 


EM3 


800 10.6 


1.5 


0.82 


3.0 


12.0 


0.65 0.33 
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700 7.8 


1.4 


0.82 


2.0 


9.0 


0.65 0.33 



hydrodynamic properties of the model ejecta, in particular, the 
density distribution in the SN envelope (Fig.[3]i. 

To estimate a measure of uncertainty in the derived phys- 
ical parameters, we investigate the sensitivity of the optimal 
model to observational values. We adopt the following rela- 
tive changes in the observational values: 20% in the bolomet- 
ric luminosity, 5% in the photospheric velocity, and 5% in the 
plateau duration. Using the auxiliary hydrodynamic models in 
the vicinity of the optimal model, we transform the adopted 
changes into changes in the pre-SN radius of ±140 Rq, the 
ejecta mass of +1 Mq, the explosion energy of +0.3 x 10^° erg, 
and the ""^Ni mass of +0.0016 Mq. Given the errors in the ob- 
servational values, these relations can be used to derive the er- 
rors in the physical parameters. The adopted relative changes in 
the observed values are close to their typical errors. We there- 
fore consider the derived changes in the physical parameters of 
SN 2005cs to represent the typical uncertainties in these values. 

3. Why non-evolutionary presupernova model? 

Although our non-evolutionary pre-SN model closely resem- 
bles the massive RSG star by the heterogeneous core-envelope 
structure, the extended radius, and the helium-core mass, it 
omits a sharp jump in density and chemical composition be- 
tween the helium core and hydrogen envelope, which is char- 
acteristic of the evolutionary model. A question may then arise: 
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Fig. 5. Optimal hydrodynamic model. Panel a): the calculated B light curve (solid line) compared with the observations of SN 
2005cs obtained by Pastorello et al. (I2006I I (open triangles) and Tsvetkov et al. (I2006I I (open circles). Panels b) and c): the same 
as panel a) but for the V and R light curves. Panel d): calculated photospheric velocity (solid line) is compared with photospheric 
velocities estimated from absorption minima of the He I 5876 A line (crosses) and the Fell 5169 A line (open circles) measured 
by Pastorello et al. (|2006] |. 



why not use evolutionary pre-SN. The answer is that the as- 
sumption of a smooth transition from the helium core to the hy- 
drogen envelope in the non-evolutionary pre-SN is dictated by 
two major facts. First, the explosion of the evolutionary model 
generally fails to reproduce the light curve of SN IIP in detail, 
as became clear after SN 1987A (cf. Wooslev [T988] l. Second, 
the Ho- profile in the SN 1987 A spectra provides clear evi- 
dence of hydrogen mixing deep down inside the helium core. 
We note that in 2D simulations the shock propagation produces 
Rayleigh-Taylor (RT) mixing between the oxygen and helium- 



core material, and between the helium-core matter and the hy- 
drogen envelope (Herant & Benz 119911 Miiller et al. 119911 
Kifonidis et al. 12003112006] !. 

In general, a self-consistent hydrodynamic model of the 
explosion of the evolutionary RSG star should be three- 
dimensional in considering both the hydrodynamic flow and ra- 
diation transfer. Unfortunately, this approach cannot presently 
be realized. We therefore accounted for 3D effects in our ID 
simulations by adopting a non-evolutionary pre-SN with den- 
sity and chemical composition jumps that had been smoothed 
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Fig. 6. Density distribution as a function of interior mass a) and 
radius b) for the evolutionary pre-SNe (Table [Tli: model EMI 
(dotted line), model EM2 (dashed line), model EM3 (thick 
solid line), and model EM4 (thin solid line). 

presumably by the RT mixing between the helium core and the 
hydrogen envelope. This approach was justified because the RT 
mixing occurs before the shock breakout and does not affect the 
light curve directly. 

Here we propose that a pre-SN has an evolutionary den- 
sity structure and chemical composition. We constructed the 
hydrostatic configuration, which reproduced the main features 
of evolutionary pre-SN models neglected before, namely the 
dense helium core with a sharp density gradient and the chem- 
ical composition jump at its boundary. We refer to this pre-SN 
model as the "evolutionary model". Four representative mod- 
els with a dense helium core and extended hydrogen envelope 
(Fig.|6]l were considered. Their parameters are listed in Table[T] 
i.e. the pre-SN radius, ejecta mass, explosion energy, total ^''Ni 
mass, helium-core mass, pre-SN mass, and surface hydrogen 
and helium abundances. The helium-core masses are typical 
of massive RSG stars (Hirschi et al. 12004 . Garcia-Berro et al. 
119971 ). The first two hydrodynamic models EMI and EM2 are 
similar to the optimal model in the basic parameters but model 
EMI differs in the initial radius. The remaining two models 
EMS and EM4 are less massive than the optimal model and 
their masses are close to the progenitor mass of SN 2005cs es- 
timated from the pre-explosion images of the galaxy M 5 1 (see 
Sect.Hl). We consider also the mixed models EMS and EM4 in 
which the helium core is mixed by the hydrogen envelope — 
in the same way as the optimal model is mixed (Fig. |2]l — in 
the inner 4.7 Mq and 2.8 Mq, respectively. 

First of all, we see that the explosion of the evolutionary 
model produces the dome-shaped light curve without a steep 
transition to the radioactive tail (Figs. |7^, b, and c). We note 
that the steep decline in luminosity at the end of the plateau 



is a characteristic of all SNe IIR The dome shape of the light 
curve for the evolutionary model is related to the almost flat 
density distribution in the hydrogen envelope (Fig.|6l). A step- 
like bump in the light curve at the transition to the radioactive 
tail is caused by the dense helium core (Fig. |6]l. We note that 
mixing in models EMS and EM4 does not remove but modifies 
the step-like feature (Figs. [TJ? and c). The fact that this bump 
is never observed in SNe IIP indicates that the density jump 
between the dense helium core and flat hydrogen envelope is 
smoothed and, consequently, the inner SN ejecta is strongly 
mixed. 

The hydrodynamic model EMI with an initial radius of 
1200 Rq, close to that of the evolutionary pre-SN (Heger 
[1998 ), predicts an unacceptably high luminosity, long plateau 
(Fig.|7^), and low photospheric velocity during the early epoch 
f < 10 days (Fig. |7}l). Model EM2 with the smaller initial 
radius, as expected, is characterized by lower luminosity, but 
remains too luminous (Fig. |7^) and has low photospheric ve- 
locity at the early epoch f < 7 days (Fig. |7}l). The reduction 
in radius would produce a reasonable luminosity at the early 
plateau stage, but this model would have too low luminosity at 
the late plateau stage. We therefore conclude that the replace- 
ment of the non-evolutionary pre-SN in the optimal model by 
the evolutionary pre-SN cannot produce a reasonable fit to the 
observational data. 

Models EMS and EM4 of lower mass and explosion en- 
ergy could reproduce the observed plateau luminosity for an 
appropriate choice of initial radius (Figs. [TJ) and c). However, 
a full description of the light curve cannot be attained for ei- 
ther the 12 Mq or 9 Mq pre-SNe. The photospheric velocity 
at the early stage in these models becomes unacceptably low 
(Fig. |7}l), which is another reason why evolutionary models 
have to be modified significantly to attain a reasonable fit to 
observational data. Remarkably, all evolutionary models fail to 
reproduce the width of the initial peak. We therefore conclude 
that there is no set of the basic parameters for the 9-18 Mq 
progenitors that can reproduce the observations of SN 2005cs 
in the frame of the ID explosion of the evolutionary pre-SNe. 

4. Progenitor mass 

To derive the main-sequence mass of the SN 2005cs progeni- 
tor, the pre-SN mass needs to be combined with the mass lost 
during the hydrogen and helium burning stages. The mass lost 
at the hydrogen burning stage is taken to be ^ 0.25 Mq, the av- 
erage value between the masses lost in the 15 Mq and 20 Mq 
evolutionary models of non-rotating stars developed by Meynet 
et al. (2003), who assumed the theoretical mass-loss rate pro- 
vided by Vink et al. ( [200TT l. 

The mass lost at the helium burning stage can be esti- 
mated from the stage duration, which is a function of the initial 
mass M and mass-loss rate M adopted for the RSG stage. We 
adopt the duration of the helium burning stage from Meynet 
et al. (i200Sj . The dependence of the mass-loss rate M on 
M was taken from Chevalier et al. (I2006I I. who used M val- 
ues from de Jager et al. (I1988I I. We also used the calibration 
M = 1 .5 X 10"^ Mq yr~' for the RSG with a main-sequence star 
M = 22 Mq, on the basis of flie SN 1999em study (Chugai et 
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Fig. 7. Bolometric light curve a), b), and c) and photospheric velocity d) of the hydrodynamic models in Table [U compared with 
the empirical data of SN 2005cs (see legends of Figs. |4]and|5]for details). Dotted line is model EMI, dashed line is model EM2, 
thick solid line is model EM3, and thin solid line is model EM4. Mixed models EM3 and EM4 are shown by the corresponding 
dotted-dashed lines in panels b) and c). 



al-EOOTll. The derived mass lost at the RSG stage was 0.6 Mq. 
The total mass lost by the stellar wind was » 0.85 Mq, and the 
main-sequence mass of the progenitor was 18.2 Mq. 

The mass of the SN 2005cs progenitor was estimated from 
archival images of the galaxy M51 taken by the Advanced 
Camera for Surveys of the Hubble Space Telescope (HST) 
and from near-infrared images acquired by the Near Infrared 
Camera and Multi-Object Spectrometer on board HST in JHK 
bands. We note that the progenitor was detected only in the I- 
band image; in other bands, only upper limits to fluxes were 
obtained. From these data, Maund et al. (I2005I I derived a pro- 



genitor mass of Li et al. (I2006I I reported a progeni- 

tor mass of 10 + 3 Mq, while Eldridge et al. ( I2007I I derived a 
progenitor mass of between 6 Mq and 8 Mq. These estimates 
therefore propose a 6 - 13 Mq range for the progenitor mass of 
SN 2005CS. 

The above mass estimates are significantly lower than our 
hydrodynamic mass. The disagreement is serious and requires 
an explanation. Our hydrodynamic model for SN IIP (Utrobin 
12004b can be checked by comparison with the independent 
model of Blinnikov et al. ( il998i ). In the case of the normal 
type IIP SN 1999em, both codes produce similar ejecta mass 
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Table 2. Hydrodynamic models for SN 1987A, SN 1999em, 
SN 2003Z, and SN 2005cs. 



SN 


Ro Mg„y E 
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(%)(Mo)(105i erg) (10-2 Mq) 
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Fig. 8. Explosion energy a) and ^''Ni mass b) versus hydrody- 
namic progenitor mass for four core-collapse SNe. 



by assuming the same distance (Baklanov et al. 120051 Utrobin 
I2007 i We explored crucial model assumptions that might min- 
imize the ejecta mass. One critical point is the degree of mixing 
between the hydrogen envelope and helium core. We found that 
the minimal mass was produced, if complete mixing occurred. 
In this case, the ejecta mass of the hydrodynamic model could 
be reduced by about 0.5 Mq. Another uncertainty was related 
to the incompleteness of the line list used in the line-opacity 
calculations. By studying this issue using the latest line list of 
Kurucz with x 6.2 x 10^ observed and predicted lines, it was 
found to provide only negligible effect, which may cause the 
mass decrease by the value of the order of 0.1 Mq. Both uncer- 
tainties implied a lower limit for the hydrodynamic progenitor 
mass as low as 17.6 Mq, which is higher than the upper limit of 
1 3 Mq recovered from the pre-explosion images of SN 2005cs. 

5. Discussion and conclusions 

The primary goal of this study was to determine parameters 
of the sub-luminous type IIP SN 2005cs by means of hy- 
drodynamic modeling. We have estimated a pre-SN mass of 
17.3 + 1 Mq, explosion energy of (4. 1 + 0.3) x 10^° erg, pre-SN 



radius of 600 + 100 Rq, and ^^Ni mass of 0.0082 ± 0.0016 Mq. 
Using conservative assumptions about the mass-loss rate at the 
hydrogen and helium burning stages, we estimate the main- 
sequence progenitor mass of 17.2 - 19.2 Mq. 

Hydrodynamic models for four SNe IIP listed in Table |2] 
are characterized by the pre-SN radius, the ejecta mass, the 
explosion energy, the total ^^Ni mass, the surface abundance 
of heavy elements, the maximum velocity of nickel, and the 
minimum velocity of the hydrogen-rich envelope. The major 
parameters of SN 2005cs — the ejecta mass, the explosion en- 
ergy, and the ^^Ni mass — are intermediate between those of 
the low-luminosity type IIP SN 2003Z (Utrobin et al. |2007] | 
and the normal type IIP SN 1999em (Utrobin l2007b in quah- 
tative agreement with their luminosities. At present, there are 
therefore four SNe IIP that have parameters determined by hy- 
drodynamic modeling. For these objects, the explosion energy 
and ^^Ni mass correlate with the progenitor mass (Fig.[8]l. This 
is consistent with the empirical relation between the explosion 
energy and -"^^-Ni mass found by Nadyozhin (2003) for normal 
SNe IIP 

Despite the uncertainties in hydrodynamic modeling, the 
disparity between the hydrodynamic mass of the SN 2005cs 
progenitor and the mass estimated from the pre-explosion im- 
ages is significant. This difference could be decreased by in- 
cluding the effects of the pre-SN light absorption in a hypo- 
thetical dusty circumstellar shell. However, this issue requires 
careful consideration, which is beyond the scope of our paper. 
We note only that this conjecture has a number of observational 
implications need to be to verified. The presence of the dense 
dusty shell around pre-SN should produce strong Na I absorp- 
tion lines in the SN IIP spectrum at the photospheric epoch. In 
the case of a normal SN IIP and normal pre-SN wind without 
a dense circumstellar shell, the predicted Na I absorptions are 
weak and probably not observable (Chugai & Utrobin l2008l l. 
In addition, the interaction of the SN ejecta with the dense cir- 
cumstellar shell should produce an outburst of radio and X-ray 
emission at an age of about ~ 10^ days. The most apparent 
effect of the light absorption in the dusty circumstellar shell 
should be a large J -K color index of the pre-SN. For instance, 
we have found that the pre-SN light absorption, required to 
allow massive progenitor implied by both the pre-explosion / 
value and upper limits in other bands for SN 2005cs, suggests a 
large color index J - K ^ 2.5 mag, compared with the intrinsic 
J -K index of 0.7 - 1 mag for typical galactic K-M supergiants 
(Elias et al. '1985i In this regard, it is noteworthy that the type 
IIP SN 2008bk in the nearby galaxy NGC 7793 with available 
JK photometry of the progenitor was found to have a moderate 
color index J - K ^ I mag (Maoz & Mannucci 2008) which 
indicates little (if any) absorption. The observational and hy- 
drodynamic studies of this supernova would be of significant 
importance in clarifying the serious and challenging problem 
of the progenitor mass of SN 2005cs and in general SNe IIP. 

The number of SNe IIP with measured hydrodynamic 
masses is too small to be able to analyze in detail and draw re- 
liable conclusions about their mass distribution. However, the 
hydrodynamic progenitor masses do appear to be systemati- 
cally higher than if SNe IIP had originated from the range of 
9-25 Mq, assuming a Salpeter initial mass function. This is 
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clearly demonstrated by the comparison of the SNe IIP mass 
distribution, calculated by assuming a Salpeter initial mass 
function in the 9-25 Mq mass range, with the mass distribution 
of four SNe IIP of known hydrodynamic mass (Fig.|9]l. Despite 
the small number of events, the significance of the difference 
between the two distributions is high: the probability that the 
four SNe occurred at random in the mass range of 15 - 25 Mq 
is only 0.01. This indicates that either the ID model of the SN 
explosion overestimates the ejecta mass, or the outcome of the 
core collapse of 9 - 15 Mq stars differs markedly from that of 
SNe IIP observed until now. 

To study the first possibility, we would require 3D radia- 
tion hydrodynamics modeling, which is not possible at present. 
While we cannot readily ascertain any 3D effect that might re- 
duce the required ejecta mass, apparent 3D effects do exist that 
could increase the hydrodynamic mass. Indeed, the RT mixing 
between the helium core and the hydrogen envelope is expected 
to produce heterogeneous ejecta consisting of helium clumps 
embedded in the hydrogen background. This structure should 
reduce unavoidably the average opacity of the hydrogen-rich 
matter As a result, the ejecta mass required to reproduce the 
observations should increase. This could counterbalance other 
possible effects that might reduce the ejecta mass. We believe 
therefore that the 3D hydrodynamic simulations are unlikely 
to reduce significantly the SN IIP progenitor masses recovered 
from the ID modeling. The situation could be eventually clari- 
fied by 3D modeling of the SN IIP explosion. 

Alternatively, the disparity between the two distributions in 
Fig.|9]is real. In this case two explanations could be invoked: 
(1)9-15 Mq stars produce faint, still undetectable SNe IIP, i.e. 
we have a selection effect; (2) core collapse of stars from this 
mass range does not produce an SN event at all. In the latter 
case, the fate of the star may be silent collapse with a neutron 
star residing inside the stellar envelope, i.e. a Thorne-Zytkow 
object (Thorne & Zvtkow ll975l) . If this is the case, the neutron 
star may eventually grow into a black hole of mass as high as 
~ 15 Mq due to the rapid (~ 10^ years) accretion of the stellar 
envelope (Bisnovatyi-Kogan & Lamzin [l9841 l. 



The second option appears, however, to be unlikely, be- 
cause it implies that, apart from double neutron stars (DNS) 
originating from 9-25 Mq stars, there should be a compara- 
ble number of binaries with a neutron star in combination with 
a black hole (NSBH binaries). Assuming that the production 
rate of DNS and NSBH binaries from 9-25 Mq stars is deter- 
mined by the random pairing of stars with the Salpeter initial 
mass function and that 9-15 Mq stars produce black holes, 
we expect the relative rate of formation of these binaries to be 
DNS : NSBH = 1 : 0.85 for 9 - 25 Mq stars. This ratio is in 
an apparent contradiction with the fact that eight DNS in the 
Galaxy are known (Ihm et al. '2006' and references therein) and 
no NSBH binary has yet been discovered. The probability of a 
random realization of this situation is only 1.5 x 10~^, i.e. suffi- 
ciently small to be able to exclude the second option that stars 
in the mass range of 9 - 15 Mq end their lifes as black holes. 

We propose that the core collapse of 9 - 15 Mq stars should 
produce a neutron star and the remaining stellar matter ejected 
as a result of a faint SN event. This picture predicts that the 
rate of faint SNe IIP should be comparable with the com- 
bined rate of normal (e.g., SN 1999em), sub-luminous (e.g., 
SN 2005cs), and low-luminosity (e.g., SN 2003Z) SNe IIP 
A detection of the extended class of faint SNe IIP, or non- 
detection at a low flux level, could verify this scenario for 
9-15 Mq stars. It is interesting that some known transient 
events, e.g. SN 1997bs (Van Dyk et al. 12000) . optical transient 
M85 OT2006-1 (Pastoreflo et al. l2007] l. and SN 2008S (Prieto 
et al. 2008 ) might belong to the proposed category of faint SNe 
related to the mass range of 9 - 15 Mq. 
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